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Abstract. We consider the distribution of cycle counts in a random regular 
graph, which is closely linked to the graph's spectral properties. We broaden 
the asymptotic regime in which the cycle counts are known to be approximately 
Poisson, and we give an explicit bound in total variation distance for the 
approximation. Using this result, we calculate limiting distributions of linear 
eigenvalue functionals for random regular graphs. 

Previous results on the distribution of cycle counts by McKay, Wormald, 
and Wysocka (2004) used the method of switchings, a combinatorial technique 
for asymptotic enumeration. Our proof uses Stein's method of exchangeable 
pairs and demonstrates an interesting connection between the two techniques. 



1. Introduction 

One of the statistics commonly studied in random matrix theory is the linear 
eigenvalue functional. If Ai, . . . , A„ are the eigenvalues of a random matrix, one 
tries to answer questions about the random variable X^ILi fi^i) some function 
/, typically as n tends to infinity. 

Suppose that G is chosen uniformly at random from the space of all simple d- 
regular graphs on n vertices, and consider its adjacency matrix. Brendan McKay 
determined the first-order behavior of its linear eigenvalue functionals, showing 
that f i^i) converged in probability to a deterministic limit |McK81| . In 

[DJPFll] , the second-order behavior of linear eigenvalue functionals was computed 
for a slightly different model of random regular graph. The motivating goal of 
this paper is to extend these second-order results to the uniform model, which we 
achieve in Theorems [12] and [T71 

We will defer further discussion of this problem and its background until Sec- 
tion |4l Until then, we discuss several combinatorial and probabilistic results that 
are interesting in their own right and that we will achieve along the way. Let Ck 
denote the number of cycles of length k in the random graph G. The distribution of 
these random variables has been studied since |Bol80[ IWor81| , where it was proven 
that (C3, . . . ,Cr) converges in law to a vector of independent Poisson random vari- 
ables as n tends to infinity, with r held fixed. As early as |McK81] . the cycle counts 
of a graph have been used to investigate properties of the graph's eigenvalues. We 
take this approach as well, converting our original problem into one of accurately 
estimating the distribution of this random vector. 
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The strongest results on the cycle counts of a random regular graph came in 
[M WWOi) ■ where the Poisson approximation was shown to hold even as d — d{n) 
and r = r(n) grow with n, so long as {d— 1)^''"^ = o(n). This is a natural boundary: 
In this asymptotic regime, all cycles in G of length r or less have disjoint edges, 
asymptotically almost surely. If {d — l)'^''^^ grows any faster, this fails. This 
led the authors in |MWW04] to speculate that the Poisson approximation failed 
beyond this threshold. Surprisingly, this is not the case. In Theorem [TTl we give a 
Poisson approximation for the cycle counts that holds so long as y/r{d — l)3/2''-i = 
o{n). We also give a quantitative bound on the accuracy of the approximation, 
which was our original motivation and is the necessary ingredient for our results on 
linear eigenvalue statistics. As a bonus, we can give in Theorem [7] a distributional 
approximation not just of the cycle counts, but of the entire process of cycles. 

The Poisson approximation in [MWW04t Theorem 1] uses a combinatorial tech- 
nique for asymptotic enumeration known as the method of switchings. We adapt 
this technique to use Stein's method of exchangeable pairs for Poisson approxi- 
mation. We discuss both methods further in the following section. As noted in 
|Wor96| ■ they have some obvious similarity, but we believe that this is the first 
time they have been connected in a rigorous way. This connection gives a novel 
construction of an exchangeable pair for use with Stein's method, and it allows the 
machinery of Stein's method to be used in some new combinatorial settings. 

In Section[2l we give some basic definitions and preliminary estimates on random 
regular graphs. Section [3] presents our Poisson approximation. The core argument 
and the most general result is Theorem [71 and our main result on cycle counts is 
Theorem 1111 In Section 21 we give the context and proofs of our results on linear 
eigenvalue functionals of random regular graphs. 

1.1. Switchings and Stein's method. The method of switchings, pioneered by 
Brendan McKay and Nicholas Wormald, has been applied to asymptotically enu- 
merate combinatorial structures that defy exact counts, including Latin rectan- 
gles |GM90| and matrices with prescribed row and column sums |McK84[ IMW031 
IGMW06] . It has seen its biggest use in analyzing regular graphs; see [KSVWOl] . 
|MWW04] ■ |KSV07j . and |BSK09] for some examples. A good summary of switch- 
ings in random regular graphs can be found in Section 2.4 of jWor99] . 

The basic idea of the method is to choose two families of objects, A and B, and 
investigate only their relative sizes. To do this, one defines a set of switchings that 
connect elements of A to elements of B. If every element of A is connected to 
roughly p objects in B, and every element in B is connected to roughly q objects 
in A, then by a double-counting argument, is approximately q/p. When 

the objects in question are elements of a probability space, this gives an estimate 
of the relative probabilities of two events. 

Stein's method (sometimes called the Stein-Chen method when used for Poisson 
approximation) is a powerful and elegant tool to compare two probability distri- 
butions. It was originally developed by Charles Stein for normal approximation; 
its first published use is |Ste72) . Louis Chen adapted the method for Poisson ap- 
proximation |Che75j . Since then. Stein, Chen, and a score of others have adapted 
Stein's method to a wide variety of circumstances. The survey paper [Rosllj gives a 
broad introduction to Stein's method, and [CDMOSj and [BHJ92| focus specifically 
on using it for Poisson approximation. 
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We will use the technique of exchangeable pairs, following the treatment in 
[CDM05) . Suppose we want to bound the distance of the law of X from the Poisson 
distribution. The technique is to introduce an auxiliary randomization to X to get 
a new random variable X' so that X and X' are exchangeable (that is, (X, X') and 
{X' , X) have the same law). If X and X' have the right relationship — specifically, if 
they behave like two steps in an immigration-death process whose stationary distri- 
bution is Poisson — then Stein's method gives an easy proof that X is approximately 
Poisson. 

Switchings and Stein's method have bumped into each other several times. For 
instance, both techniques have been used to study Latin rectangles |Ste78[ [GM90| . 
and the analysis of random contingency tables in jDS98j is similar to combinatorial 
work like [GM08j . Nevertheless, we believe that this is the first explicit connection 
between the two techniques. The essential idea is to use a random switching as the 
auxiliary randomization in constructing an exchangeable pair. 

We believe the connection between switchings and Stein's method may prove 
profitable to users of both techniques. Using Stein's method in conjunction with a 
switchings argument allows for a quantitative bound on the accuracy of the approx- 
imation. Stein's method can also be used for approximation by other distributions 
besides Poisson, and for proving concentration bounds (see |Cha07| '). On the other 
hand, Stein's method cannot prove results as sharp as |MWW04l Theorem 2], which 
gives an extremely accurate bound on the probability that a random graph has no 
cycles of length r or less. The bare-hands switching arguments used there might be 
useful to anyone who needs a particularly sharp bound on a Poisson approximation 
at a single point. 

Acknowledgments. The author gratefully acknowledges loana Dumitriu for point- 
ing out similarities between .MWW04. and [DJPPll] . and Soumik Pal and Elliot 
Paquette for their general assistance. 

2. Preliminaries 

A d-regular graph is one for which all vertices have degree exactly d. We call 
a graph simple if it has no loops (edges between a vertex and itself) or parallel 
edges. By random d-regular graph on n vertices, we mean a random graph chosen 
uniformly from the space of all simple d-regular graphs on n vertices (unless we 
specifically refer to another model). When d is odd, we always assume that n is 
even, since there are no d-regular graphs on n vertices with d and n odd. By cycle, 
we mean what is sometimes called a simple cycle: a walk on a graph starting and 
ending at the same vertex, and with no repeated edges or vertices along the way. 
For vertices u and w in a graph, we will use the notation u ^ v to denote that the 
edge uv exists. The distance between two vertices is the length of the shortest path 
between them, and the distance between two edges or sets of vertices is the shortest 
distance between a vertex in one set and a vertex in the other. 

Here and throughout, we will use ci, C2, . . . to denote absolute constants whose 
values are unimportant to us. 

Proposition 1. Let G be a random d-regular graph on n vertices, with d < n^l'^ . 
(a) Let a he a cycle of length k < n^^^^ in the complete graph Kn. Then 

P|„cq<£;<i^. 
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Figure 1. A 4-cycle and a 5-cycle, connected by a path of length 3. 



(b) Let (3 be another cycle in Kn of length j < n^^^^, and suppose that a and (3 
share f edges. Then 

co(d- 

(c) Let H be a subgraph of I'Cn consisting of a j -cycle and a k-cycle joined by path 
of length I, as in Figure[Ji Suppose that j, k, I < n^^^^ . Then 

co(d- lV+^+^ 

J>[H C G] < 'i, . 

I- — J — ^j+k+i 

Proof. These statements all follow directly from Theorem 3a in |MWW04] . □ 

3. POISSON APPROXIMATION OF CYCLE COUNTS BY StEIN'S METHOD 

3.1. Stein's method background. Recall that the main idea of Stein's method of 
exchangeable pairs is to perturb a random variable X to get a new random variable 
X' , and then to examine the relationship between the two. The basic heuristic is 
that if {X, X') is exchangeable and 

P[X'=X + 1 -, 

c 

P[X' = X-1 \X]^^, 

for some constant c, then X is approximately Poisson with mean A. (When X and 
X' are exactly Poisson with mean A and are two steps in an immigration-death 
chain whose stationary distribution is that, these equations hold exactly.) The 
following proposition gives a precise, multivariate version of this heuristic. Recall 
that the total variation distance between the laws of two random variables X and 
Y taking values in N = {0, 1,2,.. .} is given by 

dTv{X,Y) ~ sup \P[X eA]-P[YeA]\. 

ACN 

Proposition 2 ([ CDMOS) Proposition 10]). Let W = {Wi,...,Wr) be a random 
vector taking values in , and let the coordinates of Z — (Zi, . . . , Zr) be indepen- 
dent Poisson random variables with EZ/j = A^^. Let W' — {W{, . . . , W^) be defined 
on the same space as W , with iW, W) an exchangeable pair. 

For any choice of a-algebra ? with respect to which W is measurable and any 
choice of constants Ck, 

r 

dTv{W, Z) < ^a-(E|Afc - CfcP[A+ \^]\-^E\Wk-CkP[A- 

k=l 
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Vo Vi V2 V3 

AHA 

Uq Wo Ui Wi U2 W2 U3 W3 

Figure 2. The change from left to right is a forward switching, 
and from right to left is a backward switching. 

— 1/2 

with S^k — min(l, 1-4A^, ) and 

K = {K = Wk + l, W, = for k<j< r}, 

= {^k = Wk-l, W, = W; for k<j< r}. 

Remark 3. We have changed the statement of the proposition from [CDMOSj in 
two small ways: we condition our probabilities on ^F, rather than on W, and we do 
not require that EWfe = Xk (though the approximation will fail if this is far from 
true) . Neither change invalidates the proof of the proposition. 

Remark 4. There is a direct connection between switchings and a certain bare- 
hands version of Stein's method. Though this is not what we use in this paper, 
it is helpful in understanding why Stein's method and the method of switchings 
are so similar. If {X,X') is exchangeable, then as explained in [Ste92[ Section 2], 
one can directly investigate ratios of probabilities of different values of X using the 
equation 

P[X = Xi] _ P[X' =Xi\X ^X2] 
F[X = X2] ^ P[X' = X2 \ X^xiY 

This technique bears a strong resemblance to the method of switchings: if we think 
of X as some property of a random graph (for example, number of cycles) and 
X' as that property after a random switching has been applied, then this formula 
instructs us to count how many switchings change X from xi to X2 and vice versa, 
just as one does when using switchings for asymptotic enumeration. 

3.2. Counting switchings. We start by defining our switchings. Besides some 
small notational differences, the definitions will be the same as those in |MWW04] . 
To avoid repetition of the phrase "cycles of length r or less," we will refer to such 
cycles as short. 

Let G be a d-regular graph. Suppose that a — vo ■ ■ ■ Vk-i is a cycle in G, and 
let Ci — fifi+i, interpreting all indices modulo k from now on. Let e[ — WiUi+i 
for 0<i<fc— Ibe oriented edges such that neither Uj nor Wi is adjacent to Vi. 
Consider the act of deleting these 2k edges and replacing them with the edges ViUi 
and ViWi for 0<i<fc — Ito obtain a new d-regular graph G' with the cycle a 
deleted (see Figure (2). We caU this action induced given by the sequences («;), 
(wi), and (wi) a forward a-switching. We will consider forward a-switchings only 
up to cyclic rotation of indices; that is, we identify the 2k different a-switchings 
obtained by cyclically rotating all sequences Vi, Ui, and Wi. 

To go the opposite direction, suppose G contains oriented paths UiViWi for < 
i <k—l such that Vi ^ i^i+i and Wi ^ Wi+i. Consider the act of deleting all edges 



Vo Vx V2 V3 




Uo Wo Ui Wi U2 W2 U3 W3 
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UiVi and ViWi and replacing them with ViVi+i and WiWi+i for aU < i < fc — 1 
to create a new graph G' that contains the cycle a ~ vq ■ ■ ■ Vk-i- We call this a 
backwards a -switching. Again, we consider switchings only up to cyclic rotation of 
all indices. 

We call an a-switching valid if a is the only short cycle created or destroyed 
by the switching. For each valid forward a-switching taking G to G", there is a 
corresponding valid backwards a-switching taking G' to G. Let Fa and Ba be the 
number of valid forward and backwards a-switchings, respectively, on some graph 
G. Using arguments drawn from |MWWQ4l Lemma 3] , we give some estimates on 
them. 

Lemma 5. Let G be a deterministic d-regular graph on n vertices with cycle counts 
{Cfe, fc > 3}. For any short cycle a C G of length k, 



Proof. The question is, with a — vq ■ ■ ■ Vk-i and Ci — ViVi+i given, how many 
ways are there to choose e^, . . . ,e'f._i that give a valid switching? There are at 
most [njfcd'^ choices of oriented edges Cq, . . . , ej,_]^, which proves the upper bound 
dl]). For the lower bound, we demonstrate a procedure to choose these edges that 
is guaranteed to give us a valid forward a-switching. Suppose that Bq, . . . , e^_]^ 
satisfy 

(a) is not contained in any short cycle; 

(b) the distance from to is at least r; 

(c) the distance from to e'^, is at least r/2; 

(d) the distance from Wi to Ui is at least r. 

Then the switching is valid by an argument identical to the one in |MWW04] , which 
we will reproduce for convenience. By (jbl, for all i, neither Ui nor Wi is adjacent 
to Vi (or to Vi' for any i'), as required in the definition of a switching. Let G' 
be the graph obtained by applying the switching. We need to check now that the 
switching is valid; that is, the only short cycle it creates or destroys is a. 

Since a shares no edges with other short cycles, its deletion does not destroy 
any other short cycles. Condition ensures that no short cycles are destroyed 
by removing Cq, . . . , ej._]^. The switching does not create any short cycles either: 
Suppose otherwise, and let /3 be the new cycle in G". It consists of paths in GOG' , 
separated by new edges in G". Any such path in G n G" must have length at least 
r/2, because 

• if it starts and ends in a and has length less than r/2, then combining this 
path with a path in a gives an short cycle in G that overlaps with a; 

• if it starts in a and finishes inW = {uq, wq, . . . , Wfe-i, Wfc-i} and has length 
less than r/2, then combining this path with a path in a gives a path 
violating condition (|b|; 

• if it starts at some e'j and ends at e'^, then it must have length r/2 by (jg) if 
i' ^ i, and by ^ if i' — i. 



(1) 



Fa < Nfec^^ 



// a does not share an edge with another short cycle, 
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Thus f3 contains exactly one path in G DG' . The remainder of /? must be an edge 
UiVi or WiVi, impossible by (Ib|, or a path UiViWi, impossible by (|dl) . 

Now, we find the number of switchings that satisfy conditions ([l|-([d|) to get 
a lower bound on Fa- We will do this by bounding from above the number of 
switchings out of the [nj^c?'^ counted in ([ij that fail each condition (jaj)-([d|). 

• There are a total of X]^=3 edges in short cycles in G. Choosing one of 
the edges eg, ... , e'j^_^ from these and the rest arbitrarily, there are at most 
[n — l]fe-irf'^~"^/c X]j=3 "^j^j switchings that fail condition (jlj). 

• The number of edges of distance less than r from some edge is at most 
2Ej=o('^ - IV - 1 = - I)'-). At most [n - l]k-id''-^kO{{d ~ ly) 
switchings then fail condition (|b|) . 

• By a similar argument, at most [n]k-id'^^^k^O{^{d 1)''/^) switchings fail 
condition (jcj). 

• By a similar argument, at most [n]k-id''^^kO[{d — ly) switchings fail 
condition ([d|. 

Adding these up and combining O(-) terms, we find that at most 
[n - l]k-id^-^k I ^ 2jCj + 0{{d - 1)^ 

switchings out of the original [n]i;d'^ fail conditions by (jlj)-([dl), establishing □ 

For backwards switchings, we give a similar upper bound, but we only give our 
lower bound in expectation. 

Lemma 6. Let G be a random d-regular graph on n vertices, and let a he a cycle 
of length k < r in the complete graph Kn . Then 

(3) S„<(d(d-1))' 



c^kid — 1) 



and 

(4) ES„ > (d(d-l))' (^1 

Proof. The question this time is given a, how many choices of oriented paths yield 
a valid switching? For any fixed a, there are at most {d{d — l)y choices of oriented 
paths, proving ([3]). For the lower bound, let B — J2j3^P^ where /3 runs over all 
cycles of length k in the complete graph. We will first show that 

[n]k{did~-l)y / _ 4fcE;^3jQ + 0(fc(d-l)'-) \ 



(5) B> 



2k \ nd 



As in Lemma O we give conditions that ensure a valid switching. Let (5 — 
f • • • Vk-i, and suppose that the paths UiViWi in G for Q < i < k — 1 satisfy 

(a) the edges ViUi and ViWi are not contained in any short cycles; 

(b) for all 1 < j < r/2, the distance between the paths UiViWi and Ui+jVi+jWi+j is 
at least r — j + 1 . 

Any choice of edges satisfying these conditions gives a valid backwards switching: 
Condition (|b|) ensures that Vi ^ Vi+i and Wi ^ Mi+i, as required in the definition 
of a switching. Let G" be the graph obtained by applying the switching. We need 
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to check that no short cycles besides /3 are created or destroyed by the switching. 
By (jaj), none are destroyed. Suppose a short cycle /?' other than (3 is created in G". 
It consists of paths in G H G' , portions of /3, and edges WiUi+i. Any such path in 
G O G' must have length at least r/2 because 

• if it starts at Ui, Vi, or Wi and ends at Ui+j, fj+j, or WiJ^j for 1 < j < r/2, 
then ([b| implies this; 

• if it starts and ends at one of Ui, Vi, and Wi, then (jsj) implies this. 

Hence f3' must contain exactly one such path. The remainder of must either be 
an edge w^Ui+i, or a portion of /?', both of which are impossible by (|b|. 

There are [n]kd''/2k choices for /?, and at most {d{d — 1))*^ choices for Ui,Wi, 
< i < k. As before, we count how many of these potential switchings satisfy 
conditions |aj) and (|b| to get a lower bound on B. By similar arguments as in the 
proof of Lemma [5l we find that at most 

r 

2[n-l]k-i{d{d-l)f'\d-l)J2jCj 

j=3 

of the switchings violate condition (jaj), and at most [n]fe_i (d(c?— 1))'° ^0(^{d~lY~^^) 
violate condition (|b|, which proves ([5]). 

By Proposition ffil (or by |MWW04[ eq. 2.2]), 

- 2k ■ 

Applying this to ([5]) gives 

By the exchangeability of the vertex labels of G, the law of Bp is the same for 
all fc-cycles /3. It follows that EB = {[n]k/2k)'EBa, proving (g]). □ 

3.3. Applying Stein's method. Rather than prove a theorem about the vector 
of cycle counts, we will give a more general result on the process of cycles in the 
graph. Let 3 be an index set of possible cycles in Kn that the random graph G 
might contain, and for a € 3, let be an indicator on G containing a. We will show 
that the entire process (/„ , a e 3) is well approximated by a vector of independent 
Poissons, with the accuracy of the approximation depending on the size of the set 
J. We will also prove a slight variant in Proposition [TU] which achieves a better 
error bound, at the expense of considering a less general process. Our result on 
cycle counts, Theorem [TTl will follow easily from this. 

Though we have no need for these process approximations in our paper, similar 
results for the permutation model of random graph have proven useful in |JP12j . 
In any event, the machinery of Stein's method allows them to be proved with no 
extra effort. 

Theorem 7. Let G be a random d-regular graph on n vertices. For some collection 
3 of cycles in the complete graph Kn of maximum length r, we define I — {la, a € 3), 
with la = 1{G contains a}. Let Z — (Za, a £ 3) be a vector of independent Poisson 
random variables, with EZa — (d ~ where \a\ denotes the length of the 

cycle a. 
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For some absolute constant cq, for all n and d,r > 3 satisfying r < and 
d<n^/^, 



Before we give the proof, we show the resuh of applying this theorem when 3 is 
all cycles of length r or less: 

Corollary 8. Let G be a random d-regular graph on n vertices, and let 3 be the 
collection of all cycles of length r or less in the complete graph Kn . Define I and Z 
as in the previous theorem. For some absolute constant C7, for all n and d,r > 3, 

<,„,i, z, < 'jii^ini, 

n 

Proof of the corollary. If r > n^/^^ oi d > n^^^, then C'r{d — l^^^^/n > 1 for a 
sufficiently large choice of C7, and the total variation bound is trivial. Thus we can 
assume that this is not the case and apply the previous theorem: 



drvil, Z) < 



y^ C6|a|(d-l)l"l+--^ 
[n]k /cek{d- 1 



|A;+r-l 



^ 2fc V nfc+i 

^«(^)- 

The strength of Theorem [7] is that one can consider a smaller set 3 of possible 
cycles and get a tighter total variation bound. For instance, if 3 is the set of 
all cycles in Kn of length r or less containing vertex 1, then I and Z are within 
0(r{d — in total variation norm. 

Remark 9. Since the cycle counts (C3, . . . , Cr) are a functional of I, this corollary 
implies that 

dTv((C3,...,a), (Z3,...,z,)) < '^^(^-1)"'"^ 

n 

where (Z3, . . . , Z^) is a vector of independent Poisson random variables with EZ^ = 
{d — l)'^/2fc. In fact, we will give a slightly better result in Theorem [TT] 

Proof of Theorem Wc will construct an exchangeable pair by taking a step in 
a reversible Markov chain. To make this chain, define a graph © whose vertices 
consist of all d-regular graphs on n vertices. For every valid forward a-switching 
with a G 3 from a graph Go to Gi, make an undirected edge in 25 between Go 
and Gi. Place a weight of l/[rt]|o,|d'"' on each of these edges. The essential fact 
that will make our arguments work is that valid forward a-switchings from Go to 
Gi are in bijective correspondence with valid backwards a-switchings from Gi to 
Gq. Thus, we could have equivalently defined © by forming an edge for every valid 
backwards switching. 

Define the degree of a vertex in a graph with weighted edges to be the sum of 
the adjacent edge weights. Let do be the maximum degree of 25 as defined so far. 
To make 25 regular, add a weighted loop to each vertex that brings its degree up to 
do. Now, consider a random walk on 25 that moves with probability proportional 
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to the edge "weights. This random walk is a Markov chain reversible with respect 
to the uniform distribution on d-regular graphs on n vertices. Thus, if G has this 
distribution, and we obtain G' by advancing one step in the random walk, the pair 
of graphs (G, G') is exchangeable. 

Let be an indicator on G' containing the cycle a, and define I' — (/^, a E 3). 
It follows from the exchangeability of G and G' that I and I' are exchangeable, and 
we can apply Proposition [2] on this pair. Define the events A+ and A~ as in that 
proposition. By our construction. 



P[Ai I G] = 



doW|Q|d'" 

Thus by Proposition [5] with all constants set to do, 



G] 



(6) dTy(I, Z)<^E 



d\- 



EE 

ae3 



We will bound these two sums. Fix some a e 3, and let \a\ = k. By LemmalHl 



Ba 

\n\kd'' 



< 



jd-l) 
\n]k 



Thus 



E 



id -I) 



\n\k 



\n\k 



4"^ 



E 



[d-iy 



Br. 



nk 



nk 



4'' 



Applying the lower bound on E_Bq from Lemma [S] then gives 



(7) 



E 



id -I)' 



Ba 



[n\k 



[n\k 



df" 



< 



c^k{d — 1 



|fe+-r-l 



n[n\k 



In bounding the other sum, we partition our state space of random regular graphs 
into three events: 

Ai — {G does not contain a}, 

A2 ~ {G contains a, which does not share an edge with another short cycle in G}, 
A3 = {G contains a, which shares an edge with another short cycle in G}. 

On Ai, we have la ^ Fa = 0. On ^2, both bounds from Lemma [5] apply, giving us 



Fa 



d'' 



< 



On ^3, we have la — ^ and Fa = 0. In all. 



E 



\n\kd'' 



< E 



2kj:'- jCj + c4k{d-iy 



7id 



3=3 



nd 

nd 



Let d be the set of all cycles of length r or less in Kn that share no edges with a. 
On the set A2, the graph G contains no cycles outside of this set (except for a), 
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and E;=3 JCj - k + E^^g |/3| Ip- Thus 

Fa 



E 



(8) 



C4fc(d- 1)'^ 



363 



nd 
nd 



P[A2] + V[A^] 



■E/„+P[^3]. 



By Proposition nfel 

(9) 

and 

(10) 



2fc2^^ _^/k'^{d-lf 
nd " \ ra'^+i 



cikid-ir^^^^^fkid-i)''+^-' 



nd 



By Proposition [30 for any /? € 3, we have ^lalp < C2{d ~ ly^^ /n^^^ . For each 
3 < j there are at most [n]j/2j cycles in 3 of length j. Therefore 



2fc 



/36a 



J=3 



2j 



(11) 



< 



7i(i 



J=3 



lfe+1 



The last term of (|S]) is the most difficult to bound. Let % be the set of short 
cycles in Kn that share an edge with a, not including a itself. By a union bound, 



(12) 



P[^3] < ^^"^Z^- 

flex 



Now, we classify and count the cycles (3 G "JC according to the structure of a U /3. 
Suppose that /3 has length j, and consider the intersection of a and /3 (the graph 
consisting of all vertices and edges contained in both a and (3). Suppose this 
intersection graph has p components and / edges. As computed on [MWW041 
p. 5], the number of possible isomorphism types of a U /3 given p and / is at most 
{2r^)P^^ / [p — 1)!^. For each possible isomorphism type of a U /?, there are no more 
than 2kn^^P^f possible choices of f3 such that a U /3 falls into this isomorphism 
class. This is because aU /3 has j + k — p — f vertices, k of which are determined 
by a. In defining /3, the remaining j — p — f vertices can be chosen to be anything, 
and the intersection of a and f5 can be rotated around a in 2k ways, all without 
changing the isomorphism class of a U /3. In all, we have shown that the number of 
j-cycles whose overlap with a has p components and / edges is at most 



(2r3)P-i 
(p-1)!^ 



2kn^-P-f. 
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For any such choice of /3, we have Elalp < C2(d — 1)-'+'^ f /vP^'^ ^ by Proposi- 
tion [H Applying this to (fT2|). 



i=3p,/>l > 

^ (2r3)P-i 2fcc2(d- 

i=3p,/>l^'^ ' 



(13) =$:o( 



fe+1 y V Tife+i 

i=3 



Combining ([9]), (HO]), ([TT]), and (HS]), we have 



E 



Applying this and ([7]) to dl]) establishes the theorem. □ 

As mentioned in Remark [9l we can apply this theorem to give a total variation 
bound on the law of any functional of I. This bound is often less than optimal, 
since this theorem fails to exploit the factors in Proposition [21 We will take 

advantage of these factors in the following proposition, and then apply this to prove 
Theorem [11] 

Proposition 10. With the set-up of Theorem^ divide up the collection of cycles 
3 into bins !Bi, . . . , 23^. Let 

Ik = ^ ] la, Zk = y ^ Zoi, 

and let = 'EZ^ . Then 

\a\ (d- l)l"l+'-i 



k=l QSSfc 

where £_k = min (l, l-4Aj, ^^^). 

Proof. Define the exchangeable pair (G, G') as in Theorem [71 and define I[, . . . ,1'^ 
as the analogous quantities in G". Define and as in Proposition [51 noting 
that 
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By Proposition [21 



<Y,^k (E|Afc-doP[A+ I G]\ 
■(d- 1)1"! 



fc=i 



fc=i 



E 



E|4-doP[A^ I G]\ 



k=l 



E 



These summands were already bounded in expectation in Theorem [71 and applying 
these bounds proves the proposition. □ 

Theorem 11. Let G be a random d-regular graph on n vertices with cycle counts 
(Cfe, fc > 3). Let {Zk, k > 3) be independent Poisson random variables with EZfc — 
[d — 1)^ /2k. For any n>\ and r,d> 3, 

dTv{{C,,...,Cr), {Z,,....Zr)) < ^sV^id-ir/^-^ 

n 

Proof. If d > v}^"^ or r > n^/^", then C'^^id— Vf''^^'^^^ jn > 1 for a sufficiently large 
choice of cs, and the theorem holds trivially. Thus we can assume that d < n^/^ 
and r < n^/^". 

Let Afc = (d — l)''/2k. With 3k defined as the set of all cycles in Kn of length fc, 
we apply the previous proposition with bins J3 , . . . , to get 

fc(d-l)'=+'-i 

fc=3 



dTv{{G3,...,Cr), {Z3,...,Zr)) <C6 5]1.4A//' 



E 



fe=3 



o 



Vfc(d- l)'^-/2+r-l 



= o 



V^(d-l)3'V2-l 



□ 



4. Eigenvalue fluctuations of random regular graphs 



Consider a random symmetric n x n matrix Xn with eigenvalues Ai > • • • > An. 
As we mentioned in the introduction, a linear eigenvalue functional is a random 
variable of the form fi'^i) ^'^^ some function /. A common problem in ran- 

dom matrix theory is to understand the asymptotic behavior of linear eigenvalue 
functionals. Typically, one shows convergence to a deterministic limit under one 
scaling (the first-order behavior), and to a distributional limit under another scal- 
ing (the second-order behavior). The prototypical example is when is a Wigner 
matrix: the first-order behavior is given by Wigner's semicircle law (see [BSlOj for a 
modern account of Wigner's result), and for sufficiently smooth /, the fiuctuations 
from this are normal [SS981 [BY05] . 

Recently, the problem of finding the fluctuations of linear eigenvalue functionals 
was considered for random permutation matrices |BAD11) . where for sufficiently 
smooth /, the limiting distribution is not Gaussian. In fact, it is an infinitely 
divisible distribution with no Gaussian part. In [DJPPll) . the same problem was 
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5 




3 



2 



* 1 



Figure 3. The walk l-j>2-j'3-^4-j>5-J>2-j>lis non- 
backtracking, but not cyclically non-backtracking. Such walks have 
a "lollipop" shape. 

considered for the adjacency matrices of random regular graphs drawn from the 
permutation model. A random 2(i-regular graph on n vertices in this model is 
formed by choosing random permutations tti , . . . , tTcj uniformly and independently 
from the symmetric group on {1, . . . , n}, and making an edge in the graph between 
vertices i and j if 7rfc(j) = j for some k. This model of random regular graphs is 
very similar to the uniform model (see |GJKW02] for a concrete justification of this 
claim). The biggest differences are that it allows graphs to have loops and multiple 
edges, and that the expected number of fc-cycles in the limit is a{d, k)/2k, where 
a{d, k) is the number of cyclically reduced words on an alphabet of d letters and 
their inverses. This can be computed (see [DJPPll] Lemma 42]) to be 



so the limiting expectation is nearly {2d — 1)^ /2k, its value in the uniform model. 
For sufficiently smooth /, the limiting distribution of linear eigenvalue functionals 
for adjacency matrices of random graphs from the permutation model is Gaussian 
if d grows to infinity with n, and non-Gaussian if d is a fixed constant. 

Our goal is to extend this result to the uniform model of random regular graph. 
We will use Theorem [TT] to estimate the number of cyclically non-backtracking 
walks, after which most of the theorems of [DJPPH] will go through with only 
small changes to their proofs. 

If a walk on a graph begins and ends at the same vertex, we call it closed. We 
call a walk on a graph non-backtracking if it never follows an edge and immedi- 
ately follows that same edge backwards. Non-backtracking walks are also known as 
irreducible. 

Consider a closed non-backtracking walk, and suppose that its last step is any- 
thing other than the reverse of its first step (i.e., the walk does not look like the 
one given in Figure Then we call it a cyclically non-backtracking walk. These 
walks occasionally go by the name strongly irreducible. 

Let Gn be a random d-regular graph on n vertices from the uniform model, and 

(n) 

let Cf. be the number of cycles of length k in G„. We define the random variable 




{2d -1)'' ~l + 2d if fc is even, 
(2^-1)^^ + 1 if fc is odd. 
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CNBW^."'' to be the number of cyclically non-backtracking walks of length k in 
Gn- Define {cl°°\ fc > 3) to be independent Poisson random variables, with C^°°'' 
having mean Xk = (d- If /2k. It wiU be convenient to define cj°°\ C^°°\ C^"\ 



^(") 
Define 



and as zero. 



For any cycle in Gn of length j, where j divides k, we obtain 2j cyclically non- 
backtracking walks of length k by choosing a starting point and direction and then 
walking around the cycle repeatedly. In fact, if d and k are small compared to n, 
then these are likely to be the only cyclically non-backtracking walks in Gn, as we 
will prove in the course of the following theorem. 



Theorem 12. 



drv ((CNBWi"\ 1 < fc < r), (CNBW[^°' 



Proof. For any measurable function / and random variables X and Y , we have 
drv (fix), f{Y))< dTv{X,Y). It follows by Theorem [U] that 



(14) 



drv \ [J2 1 < fc < r j , (CNBW^""^ 1 < fc < r) | < 

To finish the proof, we will show that 
(15) ( H2jCj"\ 1 < fc < r j = (CNBW^"\ 1 < fc < r) 

^ 3\k ^ 

with high probability. These two vectors differ exactly when either of the following 
occur: 

Event E\. Gn contains a j-cycle and a fc-cycle with a vertex in common, with j + k < r. 
Event E2'. Gn contains a j-cycle and a fc-cycle whose distance is with / > 1 and 
j + k + 2l < r (see Figure [J). 

We have ahead done most of the work in bounding the probability of event Ei . Let 
a be some arbitrary fc-cycle. In p^ . we bounded the probability that G„ contained 
a and another cycle sharing an edge with a. With the same notation and nearly 
the same analysis (the only real change is allowing / to be zero), 

f>Q 

To bound the probability of E2, first observe that the number of subgraphs of 
Kn consisting of a j-cycle and a fc-cycle (which do not overlap) connected by a 



(00) .-^v^id-iT!^ 
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path of length I is [n]j+fe+;_i/4. By Proposition!^ each of these is contained in 
Gn with probabihty at 0({d- l)■''+''+7?^^+''+') • By a union bound, 

j+k+2l<r 

\ n 

Thus (jlSp holds with probability 1 — 0(((i — l)''/n) . If two random variables are 
equal with probability 1 — e, then the total variation distance between their laws 
is at most e. Thus the two random vectors in (|15p have total variation distance 
0{{d-iy/n). This fact and ^ prove the theorem. □ 

To relate Theorem [12] to the eigenvalues of the adjacency matrix of Gn, we define 
a set of polynomials 

ro(a;) = 1, 



d-2 



T2k{x)^2T2k[^)+j^^ forfc>l, 

T2k+i{x) ^ 2T2k+i forfc>0, 
with {r„(a;)}„gN the Chebyshev polynomials of the first kind on the interval [—1, 1]. 

Proposition 13 f [DJPPlTl Proposition 33]). Let An be the adjacency matrix of 
Gn, and let \i > ■ ■ ■ > \n be the eigenvalues of {d — 1)^^/^A„. Then 

n 

J2^k{X^) = (rf-l)-'=/2cNBWi"^ 

i=l 

By Theorem II 2[ we know the limiting distribution of X]r=i fi^i) when f{x) ~ 
Tk{x). The plan now is to extend this to a more general class of functions by 
approximating by this polynomial basis. We note the following bounds on the 
eigenvalues of uniform random regular graphs. 

Proposition 14. Let Gn be a uniform random d-regular graph on n vertices. 

(a) Suppose that d > 3 is fixed. For any e > 0, asymptotically almost surely, all 
but the highest eigenvalue of Gn is bounded by 2\/ rf — 1 + e. 

(b) Suppose that d — d{n) satisfies d = o(n^/^). Then for some constant cio, 
asymptotically almost surely, all but the highest eigenvalue of Gn is bounded by 
cioVd. 

Proof. It is well known that (jg) follows from the results in |Fri08| by various conti- 
guity results, but we cannot find an argument written down anywhere and will give 
one here. When d is even, it follows from [FriOSl Theorem 1.1] and the fact that 
for fixed d, permutation random graphs have no loops or multiple edges with prob- 
ability bounded away from zero. This implies that the eigenvalue bound holds for 
permutation random graphs conditioned to be simple, and |GJKWd2l Corollary 1.1] 
transfers the result to the uniform model. When d is odd (and n even, as it has 
to be), we apply [FriOSl Theorem 1.3], which gives the eigenvalue bound for graphs 
formed by superimposing d random perfect matchings of the n vertices. These 
are simple with probability bounded away from zero, and [Wor991 Corollary 4.17] 
transfers the result to the uniform model. 
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Fact (0 is proven in a more general context in |BFSU99l Lemma 18]. □ 

Following this definition from approximation theory, we will state the main result 
on the limiting distribution of linear eigenvalue functionals. 

Definition 15. Let p > 1, and consider the image of the circle of radius p, centered 
at the origin, under the map f{z) = — . We call this the Bernstein ellipse of 
radius p. The ellipse has foci at ±1, and the sum of the major semiaxis and the 
minor semiaxis is exactly p. 

Theorem 16. Fix d > 3, and let Gn be a random d-regular graph on n vertices 
with adjacency matrix An- Let Ai > • • • > A„ be the eigenvalues of (d — 1)~^/^A„. 

Suppose that f is a function defined on C, analytic inside a Bernstein ellipse of 
radius 2p, where p ^ (d — 1)" for some a > 3/2, and such that \f{z)\ is bounded 
inside this ellipse. Then f{x) can be expanded on [—2,2] as 



k=0 

and Y^"'^ :— X]"=i fi^i) ~ '^'^o converges in law as n —>■ oo to the infinitely divisible 
random variable 



Proof. The following facts about the Chebyshev approximation follow just as in 
Lemma 34 of jPJPPllj : 

(i) The Chebyshev series approximation for f{x) converges pointwise on the in- 
terval [2,d/^/d~^\. 

(ii) The series converges uniformly on [—2 — e, 2 + e], for some e > 0. In fact, 
defining the partial sum fm{x) = J2T=o '^k^kix), 

sup <M(d-l)-"'™, 

|a:|<2+c 

where M is a constant depending on / and d, and 3/2 < a' < a. 

(iii) The coefficients obey the bound 

lafcl <M(d-l)-"'=. 
The sum defining Yj converges almost surely, since it can be rewritten as 



oo oo 

(oo) 

i 

2 



1=1 i=i ^ ' 



and this is a sum of independent random variables, bounded in L by fact (pH)). For 
some /3 < l/a, define 

/31ogn 
log(d-l) 

4"^=E(^7^CNBWi"). 

We will use x'^^^ to approximate i^^^"'', noting that x'^^^ — /r„(Ai) — na^. By 

Theorcm[T21 the total variation distance between xl"-* and l)^'°/^afcCNBW 
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vanishes as n tends to infinity. This sum converges ahnost surely to Yf as n tends 
to infinity, so ^j-"' converges in law to Yf. By Slutsky's Theorem, we need only 



show that Y, 



in) 



f 



xj"'' converges to zero in probability. 



Fix S > 0. We need to show that 



lim P 



><5 



= 0. 



We have 



The top eigenvalue Ai is always equal to d/2\^d — 1, and by fact (P, we have the 
deterministic hmit /r„(Ai) — >■ /(Ai). Thus /(A^) — /r„(Ai) < 6/2 for all sufficiently 
large n. 

Suppose that the remaining eigenvalues are contained in [—2 — e,2 + e]. By 
fact 



(16) 



E 

1=2 



|/(A,)-/.„(A,)| <M(n-l)(d-l)-"''^" <Mn-"'^+\ 



and this tends to zero since a' (3 < 1. For sufficiently large n, this sum is thus 
bounded by S/2. We can conclude that for all large enough n, 



X 



in) 



> s 



< p 



sup |Ai| < 2 + e 

2<i<n 



and this tends to zero by Proposition [T4fel 



□ 



We can also prove that the limiting distribution of linear eigenvalue functionals is 
normal when the degree of Gn grows with n, though the conditions of the theorem 
are somewhat messier. The following theorem can be applied only when the degree 
of the graph grows more slowly than any positive power of n. This does not appear 
explicitly in the statement of the theorem, but its conditions cannot be satisfied 
otherwise. 

To remove dependence on d from our polynomial basis, define 



1 



Mx) = 2n (2 



for fc > 1. 



Theorem 17. Let Gn be a random dn-regular graph on n vertices, with dn — > 00 
as n ^ 00. Let Ai > • • • > A„ be the eigenvalues of {dn — Suppose f is 

an entire function on C, and recall cio from Proposition^!^ The function f admits 
the absolutely convergent expansion f{x) — 'Y^'^o^a'^iix) on [— cio,cio]. Denote 
the kth truncation of this series by fk := 2Zi=Q o-i^i- Let 

(3\ogn 
log(d„ - 1) 

with /3 to be specified later. Suppose that the following conditions on f hold: 
(i) For some a > 3/2 and M > 0, 

sup \f{x) - fk{x)\ < M exp{~akh{k)), 

\x\<cio 



EXCHANGEABLE PAIRS, SWITCHINGS, AND RANDOM REGULAR GRAPHS 



19 



where h is some function such that h{rn) > log(d„ — 1) for some choice of 
j3 < 1/a, for sufficiently large n. 

(n) 



lim 



= 0. 



Let fik{d) = ECNBW[.°°\ notmg that CNBW^°°^ depends on d. We define the 
following array of constants, which we will use to recenter the random variable 

mfiri) := nao + ^ , _ {nk{dn) - {dn ~ 2)nl{k is even}) 
k=l ^ " ' 

Then, as n ^ oo, the random variable 

n 

converges in law to a normal random variable with mean zero and variance <Jf — 



Proof. As n tends tends to infinity, so does r„, since assumption ^ could not be 
satisfied otlierwise. Define 



and 



noting that xj"^ — Yll=i frS^i) ~ '^/(")- (Also, note that CNBW[,°°^ depends on 

dn-) 

Let 

^(„) ^Udn- l)-'=/2(cNBWi°°^ - ECNBWi°°^) if fc < r„, 
^ I otherwise, 



and let Z3, Z^, ... be independent normals with EZ^ = and EZ^ = 2k. We will 
show that {n'^\ fc > 3) converges in law to [Zk, fc > 3) as n — > cx). Rewrite n'^^ 



as 

i<fc 

The first term converges to a centered normal with variance 2fc as n — >■ 00, by the 
normal approximation of the Poisson distribution. The random variables C^^"^ , c|°°'' , . 
are independent, so the convergence of (iV^"\ fc > 3) follows if we show that the 
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remaining terms converge to zero in probability. This holds by Chebyshev's in- 
equality, since 



Var 



It follows by the continuous mapping theorem that ^j-"'* converges to normal 
with variance af. By Theorem 1121 the total variation distance between xl"' and 



j<k j<k 

and this vanishes as dn grows 

xj-"'' approaches zero as n ^ oo, so X^"' converges in law to the same limit. 

All that remains is to show that F^"-* —X'^p'' converges to zero in probability. This 
follows exactly as in Theorem ll61 using assumptions (i) and (ii) and Proposition ll4bl 

□ 

Remark 18. The only difference between the limiting distributions of Theorems [16] 
and [T7] and those of the permutation model of random graph in |DJPP11] derive 
from the slightly different expectations of CNBW[,°°'' in the two models, and from 
the fact that CNBW^°°^ = CNBW^°°^ = in the uniform model. 
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